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ABSTRACT 

Although Lyman continuum (LyC) radiation from star-forming galaxies likely drove the reionization 
of the Universe, observations of star-forming galaxies at low redshift generally indicate low LyC escape 
fractions. However, the extreme [Oiii]/[Oii] ratios of the z = 0.1 — 0.3 Green Pea galaxies may be 
due to high escape fractions of ionizing radiation. To analyze the LyC optical depths and ionizing 
sources of these rare, compact starbursts, we compare nebular photoionization and stellar population 
models with observed emission lines in the Peas' SDSS spectra. We focus on the six most extreme 
Green Peas, the galaxies with the highest [O in] /[On] ratios and the best candidates for escaping 
ionizing radiation. The Balmer line equivalent widths and Hei A3819 emission in the extreme Peas 
support young ages of 3-5 Myr, and He n A4686 emission in five extreme Peas signals the presence of 
hard ionizing sources. Ionization by active galactic nuclei or high-mass X-ray binaries is inconsistent 
with the Peas' line ratios and ages. Although stacked spectra reveal no Wolf-Rayet (WR) features, we 
tentatively detect WR features in the SDSS spectra of three extreme Peas. Based on the Peas' ages 
and line ratios, we find that WR stars, chemically homogeneous O stars, or shocks could produce the 
observed Hen emission. If hot stars are responsible, the Peas' optical depths are ambiguous. However, 
accounting for emission from shocks lowers the inferred optical depth and suggests that the Peas may 
be optically thin. The Peas' ages likely optimize the escape of Lyman-continuum radiation; they are 
old enough for supernovae and stellar winds to reshape the interstellar medium, but young enough to 
possess large numbers of UV-luminous O or WR stars. 

Subject headings: Galaxies: evolution — Galaxies: starburst — intergalactic medium — ISM: general 
— Radiative transfer — Stars: massive 



1. INTRODUCTION 

^ . Radiative feedback from massive stars has powerful 

■ effects on both the interstellar medium (ISM) and in- 
C*") ' tergalactic medium (IGM). In the early Universe, mas- 
*/~) , sive stars likely produced most of the ionizing radia- 

■ tion responsible for r eionization (e.g.jBolton et al.ll2005t 
IBouwens et~aTJl2012t IFontanot et alJl2012f > In addition. 

, high-energy radiation from massive stars is a key com- 
^C) • poncnt of galactic energy budgets and has important 
consequences for the ISM. This radi ation may suppress 
or trigger additional star formation (Elmcgrc en fc Ladal 
Il977t iHaiman et aX1ll997l ) and may escape from Hn re- 
K> ' gions to generate the diffuse, warm ioniz ed medium (e.g., 
•jjj ; IRevnoldslll98l IQev fc Kennicuttill99l . 

The fraction of ionizing radiation that escapes H II re- 
gions and galaxies and the ionizing spectra of massive 
stars are highly uncertain, however. Near the redshift 
of reionization, the low-mass galaxies that may con- 
tribute most of the ion izing radiation are unobservable 
(|Finkelstein et al. 2012). Estimates of the galactic es- 
cape fraction necessary to sustain reionization range 
from < 13% to > 50%, depe nding on the luminosit y 
function and redshift assumed (Finkclst ein et al.l l2012). 
At low redshift, observations of star- forming galaxies 
give upper limits on the escape fraction of only a 
few percent (e.g^ iLeitherer et al.l [19951 : iHeckman et al.l 
l200lt iSiana et all 12010ft . At intermediate redshift, the 
incidence of detected Lyman continuum radiati on is 
rare (e.g., iSteidel et all l200lt IShaplev et aHl200l and 
made uncertain by p ossible low-redshift contamination 
(jVanzella et al.l 120121 ). The intrinsic ionizing spectrum 
of massive stars is likewise poorly known at all red- 



shifts. This uncertainty translates directly to uncertain- 
ties in Hn re gion and galactic escap e fractions (e.g., 
IShaplev et all I2006L iVoges et al.l 12008ft . nebular condi- 
tions, and star formation rates. 

Nebular spectra can provide information on both the 
spectral shape of t he incident UV radiation and the 
optic a l depth (e.g., IQev et al.l l2000t iGiammanco et all 
2005; IZastrow et al.l 12013ft but these observations are 
still out of reach for most high redshift galaxies. In- 
stead, many studies have focused on samples of low- 
redshift galaxies with properties similar to galaxies at 
hig h redshift. One such sam ple is the 'Green Pea' galax- 
ies (jCardamone et a l. 2009), a collection of rare emission- 
line galaxies at z = 0.1 — 0.3 identified in the Sloan Dig- 
ital Sky Survey (SDSS) through their unusually strong 
[O in] A5007 emission. In particular, due to their ex- 
tremely high [O ill] /[On] ratios, which may indicate a 
deficit of low-ionization emission, the Green Peas are 
an ideal sample to search for optically thin galaxies. 
The 80 star-forming Peas resemble high-redshift galaxies 
through their low metallicities, low extinction, high UV 
luminosities, and enormous specific star formation rates 
JCardamone et al.l2009Hlzotov et aIj201lUAmorin et al.l 
l2012aft . Hubble Space Telescope (HST) imaging of a 
handful of the Pea s reveals clumpy morp hologies and 
super star clusters (jCardamone et all 120091 ) , a common 
mode of star formation at hi gh redshift (e.g..lCowie et al.1 
[19951: [Dickinson et al.ll200l lElmegreen et al.ll2005ft . The 
Green Pea sample therefore represents the best oppor- 
tunity in the local Universe to study high redshift star- 
forming conditions and the escape of ionizing radiation. 

The precise origin of the Peas' extreme [O m]/[0 n] ra- 
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tios is uncertain. One possibility is that the high 
ratios are the result of a high ionization parameter, 
U. At a given metallicity , [O in ] /[On] increases with 
U (e .g., iKewlev fe Dopital 120021 : iMartin-Manion et al.1 
I2010D . The dimensionless ionization parameter for pho- 
toionized gas is defined as 



U 



Q 



(1) 



where Q is the hydrogen- ionizing photon emission rate, 
riH is the total hydrogen density, rp is the inner nebular 
radiu s, and c is the speed of light (|Osterbrock fe Ferlandl 
2006). As powerful starbursts, the Peas likely have com- 
pact star forming regions and a hard ionizing radiation 
field, which may act to increase their ionization param- 
eters an d explain their high [Ojn] / [O n] ratios. Alter- 
native ly, iGiammanco et ahl ()2005l ) and iPellegrini et al.1 
(2012) demonstrate that escaping ionizing radiation 
decreases the emission contribution of low-ionization 
species. This suppression would also act to increase 
the [O in] /[On] ratios and would lead to overestimates 
of U. In this paper, we evaluate the possibility that 
a low optical depth contributes to the Peas' high ob- 
served [0lil]/[0n] ratios. We specifically examine the 
Peas with the highest [O in] / [O n] ratios and consider 
whether extreme ionization parameters alone are suffi- 
cient to explain the observed ratios. 

Several previous works have investigated the ionization 
conditions in the P eas. In a dditio n to the Peas' power- 
ful [O in] emission, iHawlevi ((2012T) identifies He n A4686 
emission in nine Peas, providing further evidence for high 
ionization conditions. The presence of He n A4686 emis- 
sion in the Peas may indicate that the Peas are in the 
middle of a short-lived p hase dominated by emission from 
Wolf-Rayet (WR) stars (|Hawlevll2012l ). Additional sup- 
port for the WR sce nario comes from observations by 
lAmorm et al.l (|2012aD . who find broad, WR stellar fea- 
tures in deep spectra of three Peas. However, given the 
small sample size, it is unclear whether these observa- 
tions are characteristic of the other Peas. 

Studies of the Green Peas and similar galaxies sug- 
gest that the Peas may indeed be optically thin. In- 
termediate redshift Lya emitters (LAEs) have similar 
[O in] /[On] ratios and metallicities to the Green Peas, 
and fits to the LAEs' s pectra imply high escape fractions 
(Nakajima et al. 2012|). The Green Pea sample overlaps 
with other nearby samples of galaxies, sel ected on the ba- 
sis of high Balmer line equivalent widths (IShim fe Chary 
2012) or high UV surface brightness (Heckman et al. 
20051) . Studies of these samples also support a low op- 
tical d epth. The sam p le of Ha emitters (HAEs) stud- 
ied by IShim fe Charvl (|2012D resemble WR galaxies in 
their UV properties and inferred ionization parameters, 
although the HAEs are younger and have lower metal- 
licities than the WR galaxies. IShim fc Charvl ()2012f ) 
propose that the high ionization parameters inferred in 
the HAEs may instead indicate a low optical depth and 
a high escape fraction of ionizin g photons. Th i s sce - 
nario has also been suggested by IQverzier et al.l (|2009l ) 
for a sample of UV-bright Lyman Break Analogs (LB As), 
which includes several of the Peas. The LBAs show en- 
hanced [Om]/H/3 ratios, and ma ny of the LBAs hav e 
unusually low ratios of Ha/UV ([Overzier et al.l [2009h . 



Altho ugh a low optical dep th could explain these prop- 
erties, [OverzieFeFal] (|2009l ) conclude, that an aged star- 
burst (> 10 Myr old), dominated by supernova heating, 
provides a better explan ation for the ob s erved spectral 
properties of the sample. IHeckman et al.l (| 2 1 If) investi- 
gated the optically thin scenario by using HST observa- 
tions of Cn A1335 to infer the optical depth of several 
LBAs, including one Green Pea. Although the C n line 
in the Pea is optically thin, IHeckman et al.l ()2011[ ) point 
out that the consistency between the galaxy's Ha, UV, 
and IR luminosities makes the optically thin hypothesis 
tenuous. 

Thus, a variety of ionizing sources and optical depth 
scenarios may explain the ionization conditions in Green 
Pea galaxies. In this paper, we consider in particu- 
lar whether a low optical depth could explain the ob- 
served line ratios in the Green Peas with the highest 
[O in] / [O III ratios. We select si x extreme Green Peas 
from the ([Cardamo ne et aLl [2009) sample on the basis of 
their [Oiii]/[Oii] ratios, as these galaxies are the most 
likely to be optically thin. We then compare nebular and 
stellar population models with the observed line ratios in 
the six most extreme Peas to evaluate the galaxies' opti- 
cal depths and ionizing sources. We discuss the spectral 
line measurements and properties of the Green Peas in 
In we analyze the observed line ratios and evalu- 
ate the potential ionizing sources. We discuss the results 
and their consequences for the optical depth of the Peas 
in 2] and summarize our results in Sj5] We assume a cos- 
mology with Ho — 70 km s^ 1 Mpc -1 , Q m = 0.3, and 
Aq = 0.7 throughout. 

2. DATA AND MEASUREMENTS 

The Green Peas are a subset of extreme, compact 
emission-l ine galaxies at red sh ifts between z = . land 
z = 0.3 (|Izotov et al.M201lD . ICardamone et al.l ([20090 
identify a sample of 103 narrow-line Green Pea galax- 
ies in SDSS Data Release 7 (DR7) through color se- 
lection criteria and signal-to-noise (S/N) cuts. The 
color selection isolates galaxies whose u — r, g — r, 
r — i, and r — z colors deviate from the normal SDSS 
galaxy and quasar popula tion due to intense [O Hi] A5007 
emission in the r band (ICardamone et al~l l2~009). Us- 
ing the iBaldwin et al.l ()1981l . BPT) diagram to sep- 
arate ac tive galactic nucl e i (AG N) and star-forming 
galaxies, ICardamone et ahl (|2009f) then identify a sam- 
ple of 80 starbursting galaxies. Several of these star- 
bursts have high [Om] AA5007,4959/[On] A3727 ra- 
tios, implying unusually hig h ionization parameters (e.g., 
IMartm-Manion et al.l 120101) . The Peas' [O ill] /[On] ra- 
tios raise the possibility that the galaxies may be opti- 
cally thin to Lyman continuum radiation. A low optical 
depth would suppress the [O n] emission, consistent with 
the high observed [O in] /[On] ratios. A less-extreme 
ionization parameter in combination with a low optical 
depth could therefore be an alternative explanation for 
the observed ratios. 

We o btained spectra fo r the 80 starbursting Peas 
listed in ICardamone et all (|200l ) from the SDSS DR7 
(|Abazaiian et al.l 120091) . We then measured the fluxes 
and equivalent widths (EWs) of the galaxies' emission 
lines using the IRAFQ task splot. As discussed by 

1 IRAF is distributed by the National Optical Astronomy Ob- 
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lAmorin et all (|2012b[ ). the Peas' emission line profiles 
have broad wings, indicating the presence of high ve- 
locity gas. When necessary, we deblended the Ha and 
[Nil] lines using a Voigt profile fit, as a Gaussian fit did 
not capture the emission line wings. The [S n] AA6716 
and 6731 lines were deblended with a Gaussian fit, if 
necessary. We calculated errors in the line measurements 
according to 

Clinc = CT rms ViV, (2) 

where er rms is the average root-mean-square value per 
pixel measured at either side of the spectral line and ./V 
is the number of pixels in the line. 

We simultaneously fit the measured ratios of Ha/H/?, 
H7/H/3, and H<5/H/3 for both reddening and stellar ab- 
sorption. We adopt ed intrinsic ratios of 2.86 , 0.468, and 
0.259, respectively ([Storey fc Hummer] 1 19951) . which are 
appropriate for n H = 100 cnr 3 and T e =10,000 K. We 
then corrected the l ine fluxes for extinction using the 
iCalzetti et al.l (|2000f) reddening law. In many of the 
Peas, the [Om] AA5007,4959/A4363 ratio indicates an 
electron temperature much higher than T e =10,000 K. As 
a result, we re-fit the reddening and stellar absorption for 
these objects, using th e Balmer line ratios at T e =12,500 
K or 7W5.000 K from lStorev fc Hurrrmeil (|l995l ). To es- 
timate the appropriate temperature, we used the IRAF 
temden routine, assuming a density of 100 cm" 3 . We 
used the errors in the Ha/H/3 ratios to estimate the un- 
certainties in the extinction correction. We further as- 
sumed an uncertainty of 0.04 in the intrinsic Ha/H/3 ra- 
tio, corresponding to the change in the intrinsic ratio 
when T e is varied by 2500 K. The intrinsic ratio is less 
sensitive to density; increasing the density by an order 
of magnitude changes the intrinsic Ha/H/3 ratio by only 
0.01 or less, depending on the assumed temperature. We 
then added the reddening u ncertainty in quad r ature with 
the flux errors. As noted bv lCardamone et ail (|2009l) . the 
Peas generally have low extinction, and our average ex- 
tinction coefficient c(H/3)= 0.27. The stellar absorption 
values are also typically low, with a median absorption 
of 1.05 A, indicating that the spectra are dominated by 
a young population. 

While the Peas generally have high [O ill] /[On] ra- 
tios, two Peas (J133940.71+552740.0 and 
J094347.22+262042.5) actually have stronger [O 11] emis- 
sion than [Om]. Since we are investigating the causes 
of high [O ni]/[0 11] ratios in the Peas, we have excluded 
these objects from further analysis. 

We select the six Peas with the highest [O m]/[0 11] ra- 
tios for more detailed analysis (Figured]). These extreme 
Peas have the greatest potential for a low optical depth 
or unusual ionizing source. Furthermore, five of these ex- 
treme Peas have detectable He 11 A4686 emission, which 
improves our constraints on their ionizing sources (©. 
Line strengths, errors, and EWs for the six extreme Peas, 
the Peas with the highest [O in] /[On] ratios, are listed 
in Tabled 

The Peas have high electron temperatures, as expected 
for low-metallicity objects. After correcting for redden- 
ing and stellar absorption as discussed above, we re- 

scrvatorics, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



calculated the electron temperature. We then calculated 
the electron density from the observed [S 11] A6716/A6731 
ratio with temden and iterated until both tempera- 
ture and density converged. The revised temperatures 
are still consistent with the temperatures assumed in 
the reddening calculations. Our calculated tempera- 
tures and densities and the metallicity measurements of 
llzotov et al.l (12011 ) are li sted in Table [2] for the extreme 
Peas. Ilzotov et a" (]201lD determined abundances using 
the direct method, which calculates the oxygen abun- 
dance given the observed [O ill] and [0 11] fluxes and the 
electron temperature derived from the [O in] (A4959 + 
A5007)/A4363 ratio. Tabled also lists the SFR, derived 
from the Ha luminosity according to 



log SFR = log L(Ha) - 41.27, 



(3) 



where SFR is in units of M©yr 
-1 ~~ 



_1 and L(Ha) is 
m ergss" + (|Kennicutt fc Evans! 12012] ). Our SFRs 
are higher than th e ones found for these objects 
bv iCardamone et al l (I2009T ) . despite our use of the 
iKennicutt fc Evansl (|2012|) relat i on, w hich leads to lower 
SFRs relative to the iKennicuttl (|1998l ) relation. We find 
Ha luminosities for the extreme Peas that are typically 
a fact or of 2 greater than found by ICardamone et ail 
(200 9|) due to our different treatment of reddening. 
ICardamone et all (|2009l ) assume an intrinsic Ha/H/3 ra- 
tio appropriate for T e = 10 4 K for all the Peas, whereas 
the six extreme Peas ha ve temperatures of ~15000 K 
(Ta ble El). In addition ICardamone et al.l (|2009D use 
thelCardelli et all (119891) reddening law, instead of the 
ICalzetti et all (|2000l ) law we assume. 

As evident in Table [H the extreme Peas have an aver- 
age metallicity near 0.2 Z . At metallicities Z < 0.2 Z , 
temperature gradients in nebulae can cause the direct 
method to under- or over-estimate metallicities by 0.2- 
0.5 dex (|L6pez-Sanchez et al.l 12011 . [O 11] temperature 
estimates are necessary to correct for this effect, but un- 
fortunately, the [O 11] AA7319 and 7330 lines are dom- 
inated by noise in the SDSS spectra. This metallicity 
uncertainty should be kept in mind for the extreme Peas 
J121903, J100918, and J030321, which have metallicities 
< 0.2Z Q . 

To identify weak lines in the Peas' spectra, we created 
stacked spectra for several subsets of the Peas. We first 
corrected the spectra for redshift using the IRAF task 
dopcor. We then normalized the spectra to the mean flux 
m the 4200-4300 A range, which is free of strong spectral 
lines. As a check, we also created stacked spectra using 
the SDSS continuum-subtracted spectra and no normal- 
ization. These additional spectra had a flat continuum, 
making the hydrogen absorption lines more easily vis- 
ible. Aside from this difference, using the continuum- 
subtracted spectra does not change our results. Follow- 
ing the normalization, we re-binned the spectra to a dis- 
persion of 1 A/pixel and averaged the spectra, rejecting 
pixels that deviated by more than 5cr from the median. 
We created stacked spectra for the full sample of 78 star- 
forming Peas with [Om] A5007 >[On] A3727 and for a 
subset of 53 Peas with the strongest [O in] emission. Al- 
though our analysis focuses on the six extreme Peas, we 
created an additional stacked spectrum of the 15 highest 
[O ill] /[Oil] ratio Peas to improve the S/N. 
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TABLE 1 

Selected Emission-Line Intensities of the Extreme Green Peas 
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0.807 ± 0.008 


0.836 ± 0.077 


Hci 3819 


0.016 


± 


0. 


003 


0.018 


± 


0. 


003 
































[Nciii] 3869 


0.511 


± 


0. 


004 


0.569 


± 


0. 


049 


0.544 


± 





.017 


0. 


603 


± 


0. 


076 


0.376 ± 0.006 


0.386 ± 0.035 


Ht5 4105 


0.241 


± 


0. 


003 


0.260 


± 


0. 


022 


0.267 


± 





.011 


0. 


288 


± 


0. 


036 


0.218 ± 0.004 


0.221 ± 0.020 


H7 4340 


0.438 


± 


0. 


003 


0.461 


± 


0. 


037 


0.475 


± 





.014 


0. 


500 


± 


0. 


059 


0.426 ± 0.005 


0.423 ± 0.018 


[Om] 4363 


0.115 


± 


0. 


003 


0.120 


± 


0. 


010 


0.124 


± 





.010 


0. 


130 


± 


0. 


019 


0.115 ± 0.004 


0.115 ± 0.010 


[Fera]+Civ 4658 


0.010 


± 


0. 


003 


0.011 


± 


0. 


003 




















0.032 c ± 0.005 


0.031 c 


± 0.006 


Hen 4686 


0.007 


± 


0. 


002 


0.008 


± 


0. 


002 




















0.013 ± 





.003 


0.012 


± 


0.003 


[Ariv]+Hei 4713 


0.021 


± 


0. 


003 


0.022 


± 


0. 


003 




















0.017 ± 





.003 


0.017 


± 


0.003 


[Ariv] 4740 


0.013 


± 


0. 


003 


0.013 


± 


0. 


003 
































[Om] 4959 


2.449 


± 


0. 


010 


2.417 


± 


0. 


182 


2.236 


± 





.032 


2. 


216 


± 


0. 


241 


2.001 ± 





.013 


1.927 


± 


0.152 


[Om] 5007 


7.392 


± 


0. 


028 


7.264 


± 


0. 


545 


6.765 


± 





.090 


6. 


679 


± 


0. 


722 


6.168 ± 





.035 


5.926 


± 


0.464 


[Nl] 5200 
















































Hci 5876 


0.127 


± 


0. 


002 


0.116 


± 


0. 


008 


0.131 


± 





.010 


0. 


121 


± 


0. 


015 


0.114 ± 





.004 


0.105 


± 


0.008 


[Oi] 6300 


0.034 


± 


0. 


002 


0.030 


± 


0. 


003 




















0.016 ± 





.004 


0.015 


± 


0.004 


[Nn] 6548 


0.018 


± 


0. 


005 


0.016 


± 


0. 


004 
































Ha 6563 


3.189 


± 


0. 


013 


2.787 


± 


0. 


185 


3.193 


± 





.050 


2. 


828 


± 


0. 


272 


3.064 ± 





.020 


2.786 


± 


0.192 


[Nil] 6584 


0.069 


± 


0. 


004 


0.060 


± 


0. 


005 


0.112 


± 





.024 


0. 


099 


± 


0. 


023 


0.078 ± 





.009 


0.071 


± 


0.009 


[S11] 6716 


0.097 


± 


0. 


004 


0.084 


± 


0. 


007 


0.132 


± 





.018 


0. 


116 


± 


0. 


019 


0.063 ± 





.010 


0.056 


± 


0.010 


[S11] 6731 


0.077 


± 


0. 


004 


0.067 


± 


0. 


005 


0.079 


± 





.016 


0. 


069 


± 


0. 


015 


0.052 ± 





.010 


0.047 


± 


0.009 


[Arm] 7135 


0.058 


± 


0. 


004 


0.049 


± 


0. 


005 


0.116 


± 





.031 


0. 


099 


± 


0. 


028 


0.044 ± 





.016 


0.039 


± 


0.014 


c(H/3) 








0.21 














0.19 














0.11 






F(H/3)° 


9.06 


± 


0. 


03 


14.63 ± 0.79 


1.74 


± 





.02 


2.70 


± 


0. 


21 


6.71 ± 





.04 


8.87 


± 


0.49 


EW(abs) b A 








0.84 














0.00 














4.67 






EW(H/3)A 








228.6 














184.9 














140.5 






EW(Hq)A 








1158.0 














1086.0 












812.5 







The first 22 rows list the fluxes relative to H/3 for the designated emission lines. For each object, we list the observed fluxes (F(X)) and the fluxes 
corrected for reddening and stellar absorption (/(A)). 

a We list both the observed flux and the flux corrected for reddening and stellar absorption. Fluxes are in 10 -15 erg s _1 cm -2 . 
b Best-fit stellar absorption equivalent width. 

c Includes nebular [Fe in] A4658 emission and broad Civ emission. 
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[0 IIl]/[0 II] Distribution 




2 4 6 8 10 12 14 

[0 III] AA5007,4959/[0 II] X3727 

Fig. 1. — A histogram of the [Oin] AA5007, 4959/[Oll] A3727 ratios in the Green Peas. The six extreme Peas are shown in black. 



TABLE 2 

Properties of the Extreme Green Peas 



SDSS ID 


z a 


[Oin] AA4959, 5007/A4363 


[Sn] A6716/A6731 


T„ 


n e 


Z h 


SFR C 










(K) 


(cm- 3 ) 


(Z ) 


(Moyr" 1 ) 


J081552. 00+215623. 6 


0.1410 


69.901 ± 6.295 


1.439 d ± 0.204 


15069" 


<206 / 


0.204 


4.9 ± 0.2 


J121903. 98+152608. 5 


0.1957 


65.473 ± 5.145 


0.897 ± 0.096 


15507 


1088 


0.158 


21.0 ± 0.8 


J101157. 08+130822.0 


0.1439 


75.320 ± 5.886 


1.132 ± 0.110 


14566 


378 


0.209 


9.6 ± 0.4 


J145735. 13+223201.8 


0.1488 


80.505 ± 6.121 


1.259 ± 0.108 


14156 


171 


0.224 


13.2 ± 0.5 


J100918. 99+291621. 5 


0.2219 


68.501 ± 9.140 


1.680 d ± 0.433 


15209 e 


<189 f 


0.166 


6.0 ± 0.3 


.1030321.4-075923 


0.1650 


68.502 ± 5.778 


1.193 ± 0.301 


15204 


272 


0.148 


10.1 ± 0.4 



"From SPSS DR7 

''From llzotov et al.l 120111) . We assume a solar met allicity of 12+log(Q/H)=8 .69 JAllende Prieto et aril200lD . 

c Using the conversion from Ha luminosity given in iKennicutt fe Evand l|2012f) . 

d Tho [Sn] ratio is at the low-density limit, and we are unable to determine a density. 

e Calculatcd assuming a density of 100 cm -3 . 

■^The upper limit is calculated from the ler errors for the measured [S II] ratio. 



G 



Jaskot & Oey 



2.1. Equivalent Widths and Age Constraints 

Although several studies have proposed an 'aged star- 
burst' hypothesis to explain the high io nization condi- 
tions in the Peas and similar galaxies (e.g. JOverzier et al.l 
2009), a number of spectral features suggest younger 
ages, particularly for the highest ionization Peas. Ta- 
bles |3] and E] compare the EWs and corresponding age 
estimates for the extreme Peas. The Balmer line EWs are 
large; the full sample of 78 Peas has an average redshift- 
corrected Ha EW of 607 A, while the high-ionization 
Peas in Table H have Ha EWs -100-1000 A greater than 
this average. T o estimate the Peas' ages, we ran Star- 
burst99 models (|Leitherer et al.lll999[ ) at Z = 0.2 Z Q for 
a 10 6 M instantaneous burst and for continu ous star 
forma tion with the SFRs in Table We used a iKroupal 
(|200l initial mass function from 0.1-150 M Q , the Padova 
evolutionary tracks and mas s-loss rates with AGB stars 
( Vazquez fc Leithererl 120051) . and the iPauldrach et al.1 
(120011 ) and iHillier fc Millerl (|1998| ) stellar atmospheres. 
For the Peas' Ha EWs, the instantaneous burst model 
places upper limits on the burst ages of 4.2-5.1 Myr (Ta- 
ble [3]). The H/3 EWs suggest slightly younger ages, with 
upper limits of 3.7-4.9 Myr. These latter values may 
be more reliable due to the low continuum levels in the 
Ha spectral region and blending with the [N n] lines. We 
note, however, that both the H/3 and Ha EWs give sim- 
ilar ages. 

An additional indication of the Peas' young ages 
comes from the detection o f the Hei A3819 line. 
IGonzalez Delgado et al.1 (|1999( ) show that this weak line 
is easily covered up by stellar absorption and is only 
present for the first 3 Myr following an instantaneous 
burst of star formation. The He I A3819 line is evident in 
the individual spectra of 5 of the 78 Peas, including Peas 
J101157and J145735, which are among the extreme Peas 
(Figure [2]). Pea J121903 may also exhibit Hei A3819 
emission, but the emission is only present at a 2.6a level. 
The A3819 line shows up clearly in each of the stacked 
spectra as well (Figure (5)). 

The H Balmer lines in the stacked spectra also point 
to a young age. As with the Hei A3819 line, stcllar 
absor ption quickly dominates the h igher order Balmer 
lines ([Gonzalez Delgado et al.l[l999f ). The stacked spec- 
tra show the entire Balmer series in emission, which again 
indicates an upper limit of 3 Myr on the age of an in- 
stantaneous burst (Gonzalez Delgad o' et al.|[l999|h 

If the star formation is continuous, however, the Peas' 
ages may be older than 3 Myr. In the case of continu- 
ous star formation Starburst99 models, the Balmer line 
EWs give upper age limits of 7-60 Myr (Table 01 ■ Sim- 
ilarly, the Balmer emission in the stacked spectra gives 
an upper age limit of 20 Myr. The observed Hei A3819 
EWs indicate young ages, however, even with continu- 
ous star formation. Although Hei A3819 can exist in 
emission for the first 10 Myr of continuous star forma- 
tion, the observed Hei A3819 E Ws imply ages of 3-6 Myr 
(|Gonzalez Delgado et aTlll999D . 

Table [3] also lists the stellar mass formed in the burst. 
To estimate the size of the burst, we compare the ob- 
served Ha luminosities with the ionizing photon pro- 
duction rate, Q(H), in the Starburst99 instantaneous 
burst model. We convert Q(H) to an Ha luminos- 
ity using the Balmer line ratios, Case B emissivities, 



and Case B recombinat ion coefficients at 15,000 K from 
iStorev fe Hummeil ([19951 ). We calculate the starburst 
mass twice: once assuming an age of 3 Myr, as suggested 
by the Hei A3819 emission, and once with the instanta- 
neous burst age from the H/3 EW. Since the Ha lumi- 
nosity is roughly constant for the first 3 Myr, the burst 
mass corresponding to an age of 3 Myr should also be 
valid for younger ages. As shown in Table [3l the derived 
burst masses range from 14-171 million M^ . Usin g the 
Peas' total stellar masses from llzotov et al.l ([20111 ) and 
the 3 Myr ages, the bursts account for 2-26% of the total 
stellar mass of the Peas. If the older ages of 4-5 Myr are 
adopted, the fraction of the total mass increases. While 
the burst fractions for four of the extreme Peas remain in 
the 2-26% range, Peas J121903 and J100918 have burst 
fractions of 76% and 87% if the older ages are used. Such 
high fractions may not be realistic, however. In particu- 
lar, given the detection of He I A3819 emission in J121903, 
the younger age and lower burst fraction are preferable. 

In either case, continuous or instantaneous star for- 
mation, the Peas are evidently young and powerful star- 
bursts. We now turn to the question of the different ion- 
izing sources that may be present in these objects and 
whether they can explain the observed emission. The 
Peas themselves may span a range of ages, and different 
sources may dominate the ionization in individual galax- 
ies. 

3. HIGH IONIZATION LINES AND POSSIBLE ORIGINS 

The Green Peas have [Oiii] AA5007,4959/[On] A3727 
ratios ranging from 1.5 to as high as 13.7, which may indi- 
cate an extremely high ionization parameter and/or a low 
optical depth. Several of the Peas show other high ioniza- 
tion lines in their spectra, with strong [Ne in] A3869 emis- 
sion as well as detectable [Anv] A4712, [Anv] A4740, 
and Hen A4686. In particular, we detect Hen emission 
in 5 of the 6 most extreme (i.e., highest [O m]/[0 n] ratio) 
Peas. The presence of these weak, high ionization lines 
suggests the existence of a hard ionizing radiation field, 
such as that produced by WR stars. Interestingly, while 
the high-ionization lines show up clearly in the stacked 
spectra of the Peas (Figure [3]) , no stellar WR features 
are visible. In particular, the blue bump at 4640-4650 A, 
red bump at 5808 A, and broad He II A4686 component 
are absent. Instead, the He II emission line appears nar- 
row, suggesting a nebular origin. In the Peas' individual 
spectra, the S/N is lower, but WR features are tenta- 
tively present in three of the extreme Peas (Figure |4|) . 

Some or all of the Peas may therefore fall into a 
mysterious cl ass of galaxies with Hen emission but no 
WR features dGuseva et al.ll2000l iThuan fc Izotovll2005l: 
[Sh irazi fc Brinchmannl 120121 ) . iShirazi fc Brinchmannl 
(2012) have recently analyzed 2865 He n-emitting galax- 
ies in SDSS. They find that the fraction of He n-emitters 
without WR feature s systematically increases w ith de- 
creasing metallicity ([Shirazi fc Brinchman^ll2012H . Pos- 
sible origins for the He II emission include WR stars with 
weak stellar lines, chemically homogeneous massive stars, 
AGNs, high-ma ss X-ray binaries (HMXBs), a nd fast, ra- 
diative shocks (Shirazi fc Brinchmannl 120121 ). The pre- 
cise cause of the He n emission may vary from galaxy to 
galaxy, however. Here, we specifically evaluate these pro- 
posed ionizing sources for the six extreme Peas. We note 
that if the Peas are indeed optically thin, their low opti- 
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Fig. 2.— The top three panels show the 3740-3860 A region for Peas J121903, J101157, and J145735. J101157and J145735have detections 
of the weak He I A3819 line in emission, and He I A3819 may be present in J121903as well. The bottom panel shows the same region for 
the stacked spectrum of the 15 Peas with the highest [O III] /[Oil] ratios. 
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TABLE 3 

Ages of the Extreme Green Peas (Instantaneous Burst) 



Galaxy 




Ha 




H/3 


Hoi A3819 








EW° 


Age" 


EW a 


Age" 


EW 


Agc<= 


M B urst(3 Myr) d 


M B u rs t(H/3 EWy 




(A) 


(Myr) 


( A ) 


(Myr) 


(A) 


(Myr) 


(10 6 M o ) 


(10 6 M o ) 


.7081552 


845 


< 4.9 


196 


< 4.4 






13.5 


39.9 


J121903 


1304 


< 4.5 


206 


< 4.4 






57.9 


170.7 


J101157 


1550 


< 4.2 


293 


< 3.7 


1.91 


< 2 


26.5 


42.0 


J145735 


1008 


< 4.8 


199 


< 4.4 


1.99 


< 2 


36.4 


107.4 


J100918 


889 


< 4.9 


151 


< 4.7 






16.6 


64.4 


,1030321 


697 


< 5.1 


121 


< 4.9 






27.7 


135.6 



"Rcdshift-corrcctcd. 

b Age using the instantaneous burst Starburs t99 model. 

c Agc using the instantaneous burst mode 1 of Gonzalez Dclgado ct al. (1999). 
d From the Ha luminosity, assuming an age of 3 Myr. 

e From the Ha luminosity, using the H^J EW-dcrivcd age from Column 4. 



TABLE 4 

Ages of the Extreme Green Peas (Continuous Star Formation) 



Galaxy 




Ha 




H/3 


He I A3819 




EW 


Age" 


EW 


Age" 


EW 


Age c 




(A) 


(Myr) 


(A) 


(Myr) 


(A) 


(Myr) 


J081552 


845 


< 37.4 


196 


< 16.0 






J121903 


1304 


< 13.5 


206 


< 13.9 






J101157 


1550 


< 9.7 


293 


< 7.5 


1.91 


< 3 


J145735 


1008 


< 25.5 


199 


< 15.3 


1.99 


< 3 


J100918 


889 


< 33.4 


151 


< 31.9 






J030321 


697 


< 59.3 


121 


< 58.6 







a Rcdshift-corrcctcd. 

Age assumin g the continuous star formation Starburst99 model. 
c Agc from the Gonzalez Dclgado ct al. (1999) continuous star formation model. 




FlG. 3. — The left panel shows the 4500-4800 A region of the stacked spectrum for the 15 Peas with the highest [Olll]/[Oll] ratios. 
The right panel displays the same region for a stacked spectrum of all 78 star-forming Peas. The nebular [Fein] A4658. Hell A4686, 
[Ariv]+Hel A4713, and [Ariv] A4740 lines are labelled. No broad WR features are visible. 
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Fig. 4. — The WR blue bump and Hell emission for the three Peas which may have WR features. 



cal depths would further enhance their He n/H/3 ratios by 
decreasing the observed H/3 emission. We therefore use 
both density-bounded and radiation-bounded models to 
evaluate the possible explanations for the He II emission 
and line ratios observed in the most extreme Green Peas. 

In Table El we list the He II luminosities and the WR 
blue bump luminosities or la limits for the 6 extreme 
Peas. We calculate the la blue bump and Hell detec- 
tion limits using Equation [21 correcting for reddening, 
and assuming rest-frame line widths of 75 A and 7.07 
A, respectively. The Hen width of 7.07 A is the aver- 
age rest-frame width of the detected A4686 lines. The 
adopted blue bump width corresponds to the maximum 
exten t of the blue bump seen in th e composite WN spec- 
tra of ICrowther fe Hadfieldl (]2006f ) . As the actual extent 
of the blue bump may be less than this value (see e.g., 
Fcrnandes et al.l l2004fl . these limits are likely overesti- 
mates. In addition, we calculate the detection limits for 
the WR Si in A4565 line by assuming it has the same 
width as the Hen line. Any proposed explanation for 



the high ionization emission must be consistent with the 
ages and starburst masses of the galaxies and these con- 
straints on the WR and He n emission. 

3.1. Stellar WR Feature 

Although the narrow width of the observed He n A4686 
line seems to indicate a nebular origin, a narrow stel- 
lar line could contribute some portion of the emis- 
sion. In particular, while most WR stars have broad 
Hen emission lines, the WN9 subtype has a Hen A4686 
full width at half maximum (FWHM) of only 6 ± 
2A (jChandar et all 120031 ). Our observed Hen FWHMs 
range from ~ 3 — 5 A, slightly lower than the WN9 
FWHMs. For comparison, the SDSS spectral resolu- 
tion at these wavelengths is ~ 2 — 3 A, which sug- 
gests these lines are unresolved and not stellar in ori- 
gin. Given the low S/N of the Hen line and the uncer- 
tainty in the WN9 FWHM, however, we cannot rule out 
the WN9 stellar line hypothesis. Furthermore, observa- 
tions suggest that late-type WR stars may have weaker 
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TABLE 5 

He n and WR Emission in the Extreme Green Peas 



Galaxy 


L(Ha) 


L(Hcn A4686) 


Blue Bump L(Nni A4640) 


L(Sini A4565)° 


L(Civ A4658)° 




(10 42 crgss" 1 ) 


(10 39 ergss" 1 ) 


(10 39 crgss" 1 ) 


(10 39 ergss- 1 ) 


(10 39 ergss" 1 ) 


J081552 


0.91 


5.87 


< 4.77" 


< 1.45 




J121903 


3.91 


28.7 


< 11.2° 


< 3.41 




J101157 


1.79 


7.13 


9.79 


< 1.98 




J145735 


2.46 


6.64 


< 7.03° 


< 2.14 


9.32 


J100918 


1.12 


< 4.12" 


< 13.4° 


< 4.08 




J030321 


1.87 


8.24 


< 5.79° 


< 1.76 


20.8 



°ltr limit. 

°May include nebular [Fc III] emission. 



or narrower lines at lower m etallicities ()Conti et al.lll989t 
ICrowther fc Hadfieldll200l . 

Assuming the He n is purely stellar in origin, we 
derive the number of WN9 stars necessary to pro- 
duce the observed Hen line luminosities. We adopt 
the A468 6 luminosities for l ate- ty pe WN (WNL) stars 
given in iBrinchmann et all (120081 ): L(A4686) = 4.3 x 
10 35 ergss" 1 for Z < 0.2 Z and L(A4686) = 2.47 x 
10 36 ergss -1 for Z > 0.2 Z Q . Since the Peas have av- 
erage metallicities of ~ 0.2 Zq, we calculate the number 
of WNL stars using both luminosities and list the derived 
numbers in Tabled Brinchmann et al. derived the A4686 
luminosities from WN5-WN6 stars, which are more well 
studied than WN9 stars. However, WN5-6 stars gener- 
ally have stronger Hen em ission than WN9 stars (e.g., 
ICrowther fc Hadfiel3[2006|) . The mean He II A4686 lumi- 
nosity of WN5-WN6 stars in the Large Magellanic Cloud 
(LMC ) is a factor of 2.4 greater th an that of WN7-WN9 
stars (jCrowther fc Hadficld 120061 ). We therefore note 
that the required number of WN9 stars may be underes- 
timated. 

In addition to the blue bump emission at A4640 from 
Nih WN9 stars show Si in A4565 emission (|Guseva et al.l 
2000). Based on the blue bump and Si in luminosities or 
la detection limits (Table [5]), we calculate the number 
of WN9 stars that could be present and list these values 
in Tab le El We adopt the WN L blue bump luminosi- 
ties of IBrinchmann et al.l ([2008) , which are 4.0 x 10 35 
ergss -1 for Z < 0.2 Z^ and 6 . 3 x 1 35 ergss -1 for 
Z > 0.2 Z Q . iSchaerer fc Vaccal (fl998l ) list WNL star 
blue bump luminosities of ~ 2 x 36 ergss -1 . Adopt- 
ing these higher values would result in significantly lower 
estimates of the number of WNL stars present. For 
the Sim A4565 lu minosity, we a d opt £ (A4565) = 6.6 x 
10 35 ergss -1 from IGuseva et al.1 (|200Q[ ). We note that 
each of these line luminosities may be u ncertain by a 
factor of two (|Crowther fc Hadfieldl 12006'!. Even within 
a given WN subtype, He II and blue bump luminositie s 
spa n a wide range of valu es (jCrowther fc Hadficld 2006), 
and IGuseva et all ()2000D admit that the Si ill A4565 lu- 
minosities are poorly known. 

Assuming the Z < 0.2 Z Q values, the extreme Peas 
would need to contain 13,000-68,000 WN9 stars to ac- 
count for the observed A4686 luminosities (Table [5]). As 
in W2.ll we compute the number of O stars in the Peas by 
comparing the observed Ha luminosities to the ionizing 
photon rates in Staburst99 models. We determine how 
many 10 6 M bursts would generate the observed Ha 
luminosities and scale the number of O stars in the Star- 
burst99 models accordingly. We note that this approach 
accounts for the WR contribution to the Ha emission. 



As before, we use both the H/3 EW-derived ages and an 
age of 3 Myr. From the estimated numbers of O stars, 
the required low-metallicity WN9 population indicates 
WR/O ratios of 0.05-0.34 (TableEJ). In all but one case, 
these high WR/O ratios are impossible to match with ei- 
ther the instantaneous burst or continuous star formation 
Starburst99 models. Furthermore, while the blue bump 
limits are consistent with this number of WN9 stars at 
the 3cr level, the Sim limits do not allow such a large 
population of WN9 stars. 

A population of WN9 stars with the Z > 0.2 Z Q lumi- 
nosities could be possible. The Peas would need 2,000- 
12,000 WN9 stars, resulting in WR/O ratios of 0.01-0.06. 
The Starburst99 models show that these WR/O ratios 
occur in an instantaneous burst between 3 and 5 Myr in 
age. The la blue bump and Sim A4565 detection lim- 
its are also compatible with this number of WN9 stars. 
As WNL stars may exist on the main seque nce (e.g., 
IMartins et aLll2008L [2001 IGrafener et all2lHll) . younger 
ages are also not in conflict with the existence of a WNL 
population. 

Although the presence of sufficient WN9 stars is possi- 
ble, it is by no means certain. The feasibility of the WN9 
scenario depends strongly on the assumed emission line 
luminosities. To account for the observed He n emission, 
the WN9 stars in the Peas must have the high Hen lu- 
minosities characteristic of stars with Z > 0.2 Z Q . While 
the exact metallicity dependence of WNL line luminosi- 
ties is not yet understood, the lower line luminosities as- 
sociated with metal-poor stars are probably more appro- 
priate for the Peas. Even these latter luminosities may 
be higher than is realistic, however. As discussed above, 
the adopted WN9 He n line luminosities may be overes- 
timated by a factor of 2-3. Finally, none of the higher 
S /N stacked spectra show any evidence of WR features 
(Figure [3]) , although the blue bump emission in the low- 
metallicity case should be almost as strong as He II A4686. 
Given the significant uncert ainties in the WN9 luminos i- 
ties and line widths (e.g., ICrowther fc Hadfieldl [2006) . 
we cannot definitively rule out the WN9 scenario. How- 
ever, we conclude that stellar emission from WN9 stars 
is unlikely to be the dominant source of the observed 
He II emission. 

3.2. Ionization by Hot Stars 
3.2.1. WR Stars 

The other possibility is that the He n is nebular in ori- 
gin rather than stellar. The clear presence of other high- 
ionization nebular lines, such as [Ariv] AA4712,4740 and 
strong [Ne in] A3869, supports this scenario. If the nebu- 
lar He II emission is due to photoionization from stars, the 
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TABLE 6 

Constraints on WN9 Stars in the Green Peas 



Galaxy Af(WN9 stars)" Af(WN9 stars)" Af(WN9 stars) c WN9/0 

from He II A4686 from N III A4640 from Si III A4565 from Ho II A4686 
Z<O.2Z Z>O.2Z Z<0.2Z 9 Z>O.2Z Z < 0.2 Z a Z > 0.2 Z 



J081552 


13600 


2380 


< 11900 d 


< 7570" 


< 2200 


0.12-0.30 


0.02-0.05 


J121903 


66700 


11600 


< 28000 d 


< 17800 d 


< 5170 


0.14-0.34 


0.02-0.06 


J101157 


16600 


2890 


24500 


15500 


< 3000 


0.13-0.19 


0.02-0.03 


J145735 


15400 


2690 


< 17600 d 


< 11200 d 


< 3240 


0.05-0.13 


0.009-0.02 


J030321 


19200 


3340 


< 14500 d 


< 9190 d 


< 2670 


0.06-0.21 


0.01-0.04 



"Number of stars needed to produce the observed He II luminosity. 
"Limit on the WN9 population from N III A4640. 
c ler limits on the WN9 population from Si III A4565. 
d From the ler limit. 



stars responsible must have sufficient emission shortward 
of 228 A. Extre mely hot W R stars (e.g., iSmith et all 
2002) or O stars (|Kudritzkill20"02T) may have spectral en- 
ergy distributions hard enough to photoionize Hen. 

To test this hypothesis, we create a grid of photoioniza- 
tion models with CLOUDY version 10.00 (jFerland et al.1 
1998) for nebul ae around s t ars of different spectral types. 
We use the ISmith et all (|2002fl model grids at Z = 
0.2 Z for main sequence O stars of spectral types 03- 
06 and for WN and WC stars with core temperatures 
of 45,000 K-120,000 K; the WR core temperatures corre- 
spond to a stellar optical depth of 10. The grids use WM- 
Basic stellar atmosphere models ( |Pauldrach et al.ll2001f) 
for th e O stars and CMFGEN models (jHillier fe Milled 
[l99cl for the WR stars. 

Since the ionization parameter is set by a combina- 
tion of nebular geometry, density, and input radiation 
field (Equation [T]), we vary these parameters in our mod- 
els. Our input stellar spectral energy distributions are 
scaled to a constant hydrogen-ionizing photon rate of 
Q(H) = 10 49 - 5 photons s " 1 , which is the Q (H) value 
of an 03 V star at 0.2 Z Q ([Smith et alj [20021. We have 
run additional models with Q(H) = 10 50 photons s -1 ; 
this change does not affect our results. For each spec- 
tral type, we vary the inner radius (rj n )of the nebula 
from 0.5 pc to 8 pc and the hydrogen density from 0.01 
to 1000 cm -3 . These inner radii correspond to a wide 
range of morphologies; r in ranges from 0.2% to 99.9% of 
the outer radius of the ionized gas (r out ). Morpholo- 
gies with r m /r out < 5% only appear in models with 
njj < 10 cm~ 3 . Extremely thin shell morphologies, with 
r in/ r out > 95%, appear in some of the high density mod- 
els and are marked with open circles in the figures. 

The gas phase metallicities match the stellar metallici- 
ties in each mod el run. We use element abundances from 
lMcGaughl(|1991f) . interpolating between Z — 0.15 Z Q and 
Z = 0.30 Z . For the abundances of all other ele- 
ments, we use the CLOUDY gas-phase ISM abund ances 
(jCowie fe Songailal 119861 : iSavage fe Sembachl 119961 ) and 
scale the ab undances o f the a -elements to the O abun- 
dance from iMcGaughl (|199lD . We set the turbulence 
to 4 kms -1 and the grain abundances to 20% of the 
CLOUDY ISM value (IMathis et alJ[l97l . Changing the 
turbulence to 40 kms -1 does not affect our results. 

We analyze the resulting line emission from both 
radiation-bounded and density-bounded scenarios by 
examining the cumulative line fluxes at different radii 
within the nebulae. The fluxes at low radii correspond 
to optically thin scenarios, since the gas has not fully 
absorbed all the ionizing photons at that radius. To 



evaluate the models, we compare the predicted and ob- 
served line ratios for [O n]A3727/H/3, [Ne m]A3869/H/3, 
HeiiA4686/H/3, [Ariv]A4740/[Arm]A7135, 
[Oiii]AA5007,4959/[Oii]A3727, and HeiA5876/H/3. 
We focus on [Ar iv] A4740 rather than the stronger 
[Ariv] A4711 line due to contamination from Hei A4713. 
While [O in] /[On] and [Ar iv]/[Ar in] are particularly 
sensitive to the ionization parameter or optical depth, 
Hell/H/3 sets strong constraints on the temperature of 
the ionizing source. As shown in Figure [5j even though 
both WR and O star models can match the other line 
ratios, only the hottest, early- type WR stars, with core 
temperatures greater than 70,000 K, can produce the 
observed He n/H/3 ratios. 

Usin g the He + -ionizing photon rates from ISmith et ail 
(2002) for the hottest WN and WC stars, we calculate 
the number of extreme early-type WR stars necessary 
to account for the observed Hell luminosities. We use 
the helium emissivities and r ecombination coefficients 
from iStorev fe Hummerl ([1995) for n e = 100 cm and 
T e = 15, 000 K. Tables [7] and [8] list the numbers of early- 
type WN (WNE) and WC (WCE) stars. The resulting 
WR/O ratios, calculated as in §3.1| are a few percent, 
which matches the expected value for a 3-4 Myr old in- 
stantaneous burst. However, if the star formation is con- 
tinuous, the Starburst99 models show that the WR/O 
ratios should never increase above 0.016 due to the con- 
tinued production of new O stars. 

While the WR/O ratios are reasonable, at least for an 
instantaneous burst, we also consider whether the con- 
straints on WR stellar features are consistent with this 
number of WR stars. The blue bump emission of WNE 
stars comes from Nv A4603-4620, and the WNE stars 
have weaker blue bump l uminosities than the WNL stars 
(jBrinchmann et ai1l2008l ). Thus, the required population 
of 2000-9000 WNE stars could easily go unnoticed in the 
Peas' spectra (see Table . WNE stars do have strong 
and broad He n A4686 emission, however. Adopting the 
Z < 0.2Z w and Z > 0.2 Z ^ WNE A4686 luminosities of 
iBrinchmann et aLI (|2008f ). we find that the broad com- 
ponent should be 6%-30% the strength of the nebular 
component, respectively. Although we do not see any 
broad He n component in the stacked spectra, we cannot 
rule out the presence of WNE stars, particularly if the 
lower metallicity luminosities are valid. 

The WCE stars should also have broad Hen emis- 
sion, blended with Cm A 4650 and Civ A4658 (e.g., 
iCrowther fc Hadfieldl l200ll . From the WCE A4686 and 
A4650 luminosities of ICrowther fc Hadfieldl (|2006l) . the 
broad Hen component should be 89% as strong as the 
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Fig. 5. — Radiation-bounded (optically thick) CLOUDY photoionization model grids for single-star Hn regions at densities of rtjj = 100 
cm -3 (top) and rtjj = 1000 cm -3 (bottom). The black crosses show the positions of the five high-ionization Peas with Hell detections and 
the la errors. The red grid corresponds to the WN star photoionization models, green corresponds to WC stars, and purple corresponds to 
main sequence O stars. Solid lines indicate constant stellar temperature; the O star effective temperatures and WR star core temperatures 
are labelled. Dotted lines indicate constant values of the inner nebular radius, from 0.5 pc to 8 pc. Thin shell models, where r- ln /r ou t > 95%, 
are marked with open circles. Only the hottest WR stars match the observed line ratios. 
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TABLE 7 

Constraints on WNE Stars in the Green Peas 

Galaxy 7V(WNE stars)" iV(WNE stars) 6 WNE/O 

from He II A4686 from N V AA4603-20 from He II A4686 







Z < 0.2Z o 


Z > 0.2Z© 




J 081552 


2100 


< 200000 


< 30000 


0.02-0.05 


J121903 


10000 


< 470000 


< 70000 


0.02-0.05 


J101157 


2600 


< 270000 


< 41000 


0.02-0.03 


J145735 


2400 


< 290000 


< 44000 


0.008-0.02 


J030321 


3000 


< 240000 


< 36000 


0.009-0.03 



"Number of stars needed to produce the observed He II luminosity. 

6 lrx limit on the WNE population from N V AA4603-4620, using the WNE star N V luminosities of IBrinchmann et all 1 1200811 . 



TABLE 8 

Constraints on WCE Stars in the Green Peas 



Galaxy 


JV(WCE stars)" 
from He II A4686 


JV(WCE stars) 6 
from C HI+C iv AA4650-4658 


WCE/O 
from Hen A4686 


J081552 


1100 




0.009-0.02 


J121903 


5200 




0.01-0.03 


J101157 


1300 


1200 


0.01 


J145735 


1200 


1900 


0.004-0.01 


,1030321 


1500 


1400 


0.004-0.02 



"Number of stars needed to produce the observed Hell luminosity. 

^Inferred WCE population from Cili+Crv AA4650-4658, using the WCE star Cm and Civ luminosities of lCrowther &c HadfielcH lf2006l) . The 
Clll+Civ emission may be contaminated by nebular [Fclll]. 
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nebular emission. The lack of this broad emission in the 
stacked spectra (Figure [3]) indicates that thousands of 
WCE stars are unlikely to be present in most Peas. One 
Pea that might plausibly have WCE stars is J030321, 
which appears to have broad AA4650-4670 emission (Fig- 
ure |4j . The strength of this emission is consistent with 
the number of WCE stars required to produce the ob- 
served He II. As with the WNE stars, Starburst99 models 
can match the inferred WC/O ratios with an instanta- 
neous burst but not with continuous star formation. 

The above calculations show that photoionization by 
early-type WR stars is possible. In reality, the nebu- 
lar gas in the Peas is likely photoionized by stars of 
many different spectral types. Consequently, we ran 
CLOUDY photoionization models with the 10 6 M Q Star- 
burst99 instantaneous burst model ($2jJ as the input 
ionizing spectrum. The nebular model parameters are 
the same as described for the single-star models, except 
for the inner radius, which we varied from 20-81 pc. 
As before, these models occasionally produce thin shell 
nebulae, where r in /r out > 95%. We denote radiation- 
bounded, thin shell models by open circles in the fig- 
ures. The other optically thick models have n n /r out be- 
tween 8.6 and 95%. The thin shell models may be valid 
if dense filaments dominate the Peas' nebular emission. 
In addition, unlike the single-star CLOUDY models, the 
hydrogen-ionizing photon rate, Q(H), varies as a function 
of the age of the starburst. 

We plot the CLOUDY and Starburst99 results for sev- 
eral emission line ratios in Figures [6]|9l Ages of ~ 4- 
5 Myr, when the burst contains a high fraction of WR 
stars, do the best job of fitting the observed He n/H/3 ra- 
tios. At all ages, the models have difficulty matching 
the observed [Oii]/H/3 ratios. Part of this discrepancy 
arises from the modeled nebular morphology. The mod- 
els treat the ionized g shell surrounding an inner 
cavity. While this is an appropriate model for a single su- 
perbubble, the ionized gas in the Peas is likely much less 
homogeneous. By adding a filling factor and treating the 
ionized gas as dense clumps in a lower density medium, 
we obtain an improved fit to the [O n]/H/3 emission (Fig- 
ures [8] and [9j see also Zastrow et al. 2013). We adopt 
filling factors of 1, 0.1, and 0.01, where the filling factor 
indicates the volume of the gas occupied by clumps of 
a given density. The remaining volume is treated as a 
vacuum. The filling factor models are still offset to lower 
[O m] /[On] ratios than we observe, which may indicate 
a need for higher temperature sources. As we discuss 
in the next section, unusually hot massive stars could 
explain this discrepancy. 

The CLOUDY models support both optically thick and 
optically thin scenarios, depending on the nebular pa- 
rameters. In general, the models with higher densities 
and larger inner radii imply lower optical depths. For 
densities n < 1000 cm -3 , the strength of the [O ll]/H/3 ra- 
tio implies a high optical depth, while the combination of 
[Nemj/H/3 and [Ar iv]/[Ar III] ratios suggests a low op- 
tical depth. For the n — 1000 cm~ 3 models, applicable 
to Pea J121903, the [On]/H/3 ratio supports a low op- 
tical depth if the inner radius is larger than 50 pc or if 
the filling factor is 0.01. The other line ratios generally 
favor low optical depths for n = 1000 cm~ 3 and an inner 
radius larger than 30 pc and high optical depths in the 



other models. We discuss the optical depth further in <gj 
Photoionization by WR stars is a viable explanation 
for the Peas' line ratios, as long as the starbursts are old 
enough to have a sizeable WR population. The clear de- 
tection of the He I A3819 line in Peas J121903, J101157, 
and J145735 (see §2.1j) calls this interpretation into ques- 
tion. In addition, the WR hypothesis depends strongly 
on our choice of WR line luminosities. While hot WR 
stars may account for part of the Hell emission in the 
high-ionization Peas, additional causes of the emission 
may be required for the youngest Peas. 

3.2.2. Homogeneously-Evolved O Stars 

If WR stars are not present, hot, low-metallicity O 
stars could perhaps be the source of the hard ioniz- 
ing photons. In the past few decades, studies have in- 
creasingly recognized that massive, rotating stars at low 
metallicity can achieve much higher effective tempera- 
tures than ordinary O stars. At low metallicity, weaker 
stellar winds make it diffi cult for stars to lose angular 
momentum (jLangerl fl998|) . and as a result, stars can 
maintain high rotational velocities. The associated ro- 
tational mixing may then result in chemicall y homo- 
geneous or quasi-homo geneous evolution (e.g., iMaederl 
119871 lYoon et all [2006). Chemically homogeneous O 
stars may have effective temperatures much higher than 
expecte d for their mass. For instance, iBrott et al.l 
(|2011al ) show that, at the metallicity of the LMC, stars 
more massive than 20 M Q may reach temperatures as 
high as 60,000-70,000 K. The highest mass stars are the 
most likely to evolv e homogeneously a nd reach such high 
temperatures (e.g., IBrott et al.ll2011al ). These stars also 
evolve more quickly than lower mass O star s and could 
be pres ent in young starbursts. For instance. IBrott et al.l 
(|2011a|) show that 60 M stars at Small Magellanic 
Cloud (SMC) metallicity can reach 60,000 K within 3.4 
Myr. These high temperature O stars may be plentiful 
in low-metallicity starbursts and could provide the nec- 
essary high-energy photons for the Peas' Hen emission. 

As a test, we co nsider the 60,000 K O star models 
of IKudritzkil (pOOl . We ran a single-star Hn region 
CLOUDY model as described in pXll for the T cff = 
60,000 K, log(L/L ) = 7.03 luminosity O star model. 
Figure [TU] shows that a 60,000 K O star continues the 
trend of increasing Hen/H/3 with effective temperature. 
The 60,000 K O star model still under-predicts the ob- 
served Heii/H/3 ratio. However, O stars with tempera- 
tures above 60,000 K could result in higher Heii/H/3 ra- 
tios. 

To determine the number of hot O stars needed to gen- 
erate the observe d Hell luminositie s, we use the ionizing 
photon rates from[kudritzki| ()2002h for the log(L/ L m) = 
7.03, 60,000 K star (see Figure 11 of lKudritzkJ^OOl and 
emissivities and reco mbination rates for T, = 15,000 K 
and n H = 100 cm" 3 (jStorev fc Hummed fl995h . Tabled 
lists the calculated numbers of hot O stars. The number 
of 60,000 K O stars required to explain the Peas' emission 
represent 0.4-3% of their total O star populations. For 
comparison, at ages less than 3 Myr, the instantaneous 
burst Starburst99 model shows that 8-18% of O stars 
have spectral types earlier than 04. Therefore, if the 
earliest spectral types evolve homogeneously, this pop- 
ulation could account for the high ionization lines ob- 
served in the Peas. Nevertheless, we caution that the ex- 
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Fig. 6.— CLOUDY grids with a 10 6 M ionizing cluster and a nebular density of njj = 100 cm . The high-ionization Peas are marked 
with black crosses indicating the la errors for the line ratios. The black grids have an inner radius of 20 pc, while the red grids have an inner 
radius of 81 pc. Solid lines indicate the age of the cluster and correspond to ages of 0, 3, 4, and 5 Myr. The arrow labeled 'DB' indicates 
the direction of decreasing optical depth (increasing density-bounding). Constant values of the optical depth are indicated by the dashed 
lines in each grid. The dashed lines show the emission as the nebular size is decreased from its full ionized extent to 20% of its original 
thickness, corresponding to a decrease in the optical depth. The arrow labeled 'Red' indicates the reddening vector; the plotted points 
would move parallel to this vector as the reddening correction is increased. Open circles mark thin shell models, where r; n /r out > 95% and 
''out is the outer radius of nebula for the radiation-bounded case. 



istence of chemically homogeneous massive stars has not 
yet been confirmed, even in local, low-metallicity galax- 
ies (e.g. JBrott et al.ll2011bl ). More research is necessary 
to establish whether homogeneous stars are present in 
the Peas. 

3.3. Active Galactic Nuclei 

Although the BPT diagram indicates the Peas are 
star-forming, many of th e galaxies lie close to the 
maximum starburst line of IKewlev et al.l (|2001|) . which 
marks a theoretical separation between starbursts and 
AGN. Given their proximity to this line, the Peas could 
have so me AGN contribution to the ir emission (Fig- 
ure EJ). [Sh irazi fe Brinchman n (2012) show that a plot 
of Hen A4686/H/3 versus [Nil] A6584/Ha results in a 
much cleaner separation between AGN and star-forming 
galaxi es. Figure [TT] shows the IShirazi &: Brinchmann! 
(|2012f) diagnostic for the Green Peas with Hen A4686 
detections. The line indicates the region where AGN con- 
tribute 10% of the observed Hen flux. All the observed 
He II fluxes lie below this line, solidly in the star-forming 
regime, confirming that AGN activity contributes negli- 
gibly to the high ionization emission in the Peas. 

3.4. High-Mass X-ray Binaries 

Following a starburst event, HMXBs could be a signif- 
icant source of hard photons. In addition, HMXB lumi- 



nosities may be higher at the low metallicities character- 
istic of the Peas. Given the increase of X-ray heating with 
metallicity, black holes may be able to accrete more gas 
at low metallicity (jThuan et al.l r2u04V Alternatively, the 
weaker winds and lower mass loss rates of low metallic- 
ity sta rs may lead to the formation of highe r mass black 
holes (|Thuan fc Izotovi [200l iFragos et all [2012t ). In 
thr ee low-metallicity bl ue compact dwarf (BCD) galax- 
ies, IKaaret et al.l (l2C)H find X -ray luminosities an order 
of magnitude higher than expected from the SFR. The 
metallicities of these sources are lower than the Green 
Peas, however, with Z < 0.07 Z Q ; HMXB luminosities at 
0.2 Zq might not show such a dramatic increase. 

While an HMXB scenario is appealing, it faces a num- 
ber of problems when applied to the Peas. The Peas' 
young ages may be at odds with the HMXB interpreta- 
tion, as the formation of HMXBs requires both the for- 
mation of a compact object and additional time for th e 
evolution of the companion star (IGhosh fc Whitdl200"l . 
The first HMXBs appear 4 Myr after a burst of star for- 
mation (|Linden et al.l [20 10). which is inconsistent with 
the upper age limits for at least one of the P eas. In ad- 
dition, the simulations of iLinden et al.l (|2010t ) show that 
HMXB numbers and luminosities have a complex depen- 
dence on metallicity. In fact, at 5 Myr after a burst, high- 
metallicity galaxies have a greater number of luminous 
HMXBs than low-metallicity galaxies. 
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TABLE 9 
Hot O Stars in the Green Peas 



Galaxy 


iV(60,000 K O stars) 
from He II A4686 


60,000 K O Stars / Total O Stars 


J081552 


1200 


0.01-0.03 


J121903 


5800 


0.01-0.03 


J101157 


1500 


0.01-0.02 


J145735 


1400 


0.004-0.01 


J030321 


1700 


0.005-0.02 
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Fig. 11. — The t op panel shows the Green Peas' positions on the BPT diagram (Baldwin ct al. 1981). The solid line indicates the 
IKewlev ct al. (2001) extreme starburst line. The Peas lie close to the dividing line between starbursts and AGN. R epresentative la errors 
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At a burst age of 4-5 Myr and Z = 0.2 Z^. lLinden et al.l 
• 20 10; predict the formation of 5 HMXBs with X-ray lu- 
minosities greater than 10 36 ergs s _1 per 1O 6 M0 burst. 
Four of these HMXBs would qualify as ultra-luminous 
X-ray sources (ULXs) with X-ray luminosities above 
10 39 ergss _1 . The ULX population peaks at this age, 
while the HMXB populat ion reaches its peak a few Myr 
later ([Linden et al.l |2010| ). Scaling these numbers to 
the calculated burst masses of the Peas, we would ex- 
pect 68-854 lumi nous HMXBs and 5 4-683 ULXs at an 
age of 4-5 Myr. IKaaret etHl (|2004f ) report the detec- 
tion of strong He n A4686 emission, with a luminosity of 
2.7 x 10 36 ergss~ 1 , around a ULX nebul a in Holmberg II; 
the g alaxy has a metallicity near 0.2 (|Moustakas et al.l 
2010). Adopting this luminosity for the Hell emission 
produced by a ULX, the extreme Peas would need 2100- 
11,000 ULXs to account for their He n luminosities. Even 
if the Peas are old enough to have an HMXB population, 
this number of ULXs is mo re than an ord e r of m agnitude 
larger than expected. The iGrimm et al.l (|2003f) relation 
between the SFR and the number of luminous HMXBs 
results in an even worse discrepancy, with only 20-89 
HMXBs expected in the Peas. We conclude that HMXBs 
are an unlikely source of the observed He II emission. 

3.5. Shocks 



iGuseva et all (|2000f ) and iThuan fc Izotovl ((20051 ) pro- 
pose fast, radiative shocks as a means of generating 
Hen emi ssion in galaxies withou t WR stars. In a study 
of BCDs, IThuan fc Izotovl (pOOl observe emission from 
Ne 4+ and Fe 4+ , which have ionization potentials even 
higher than He + . Radiative shocks can naturally ac- 
count for the existence of these highly ionized emission 
lines. In addition, the [Olli] AA4959,5007 and R/3 lines 
in the BC Ds exhibit broad wings indicative of high ve- 
locity gas (|Thuan fe Izo tov 2005). Observations of the 
BCD SBS 0335-052E show that the Hen and [Ne v] emis- 
sion originate i n a separate region from the [Ar iv] and 
[Oni] emission (jlzotov et al.ll200q) . The latter emission 
lines are likely photoionized by the youngest clusters in 
the galaxy. The Hen emission, on the other hand, is 
associated with older clusters, which supports an inter- 
pretation of supern o va remnant (SNR ) shoc k ionization 
(jlzotov et al1[2006h . IThuan fe Izotovl ((2005) argue that 
fast, radiative shocks are most effective in dense, low- 
metallicity, super star clusters, conditions which likely 
characterize the star-forming environments of the Green 
Peas. 

Supernov ae, galactic outflows, or merge rs could lead to 
fast shocks (|Shirazi fe Brinchmann| [2012;). and the Peas 
show evidence for these processes. Broad emission-line 
wings and multiple velocity components indicate high- 
velocit y inflows and outflows in the Peas (jAmorm et al.1 
I2012bl) . The Peas' offset in the mass-metallicity rela- 
tion likewise suggests the presence of low-metallicity gas 
inflow s or metallicity- dep endent, supernova-driven out- 
flows (|Amorin et al.ll20"Tol) . Merger-induced shocks may 
also be present, as shown by the disturbed morphologies 
and p resence of companion galaxies in HST images of the 
Peas (jOverzier et al.ll2008t ICardamone et al.ll2009ft . 

We detect emission from [Oi] A6300 and [Ni] A5200, 
which are characteristic of shocks. Weak [Ni] emis- 
sion is visible in Pea J101157's spectrum as well 



as in the stacked spectra. We do not detect the 
Ne v] or [Fev] lines discussed by IThuan fe Izotovl 
2005), but their data show that these lines are sub- 
stantially weaker than Hen. Another non-detected 
shock line is [Mgi] A 4 571. The Mappings III shock 
models of lAllen et al.1 (j2008| ) show that low ratios of 
[Mgi]/Hen can only be obtained for high magnetic 
field strengths, >10fiG. Magnetic field strengths of 20- 
50^G have been observed in starbur st galaxies (e.g., 
iKlein et all 119881 : iChvzv fc Beckll20"ol . Radio observa- 
tions indicate that the Gr een Peas have similar ma gnetic 
field strengths of > 30/iG (Chak raborti et aLll2012f ). con- 
sistent with the non-detection of [Mgi]. Figure [T2l shows 
the detected and non-detected shock lines in the high- 
ionization Pea stacked spectrum. 

Fro m the Mappings III shock model grids ([Allen et al.1 
l200l ). the observed Heii/[0 i] ratios in the extreme 
Peas support shock velocities between ~ 200 and 
~ 700 kins" 1 . The observed Ha full-width-at-zero- 
intensity values show that the Peas do have gas mov- 
ing at these velocities. In the extreme Peas, these broad 
wings range from ~ ±400-800 kms" 1 , simi lar to the ve- 
locities observed bv IThuan fc Izotovl ()2005l ). We do not 
detect any significant velocity shift in the centers of the 
Hen or [0 i] lines, and their FWHM values are similar to 
the Ha FWHM. Higher resolution observations are nec- 
essary to determine whether the Hen emission has any 
unusual velocity components. 

Both shocks and photoionization may contribute to the 
observed nebular emission. To analyze the ionizing spec- 
trum and the optical depth, therefore, we need to sub- 
tract the shock contribution. We use the lAllen et al.l 
(2008 ;) Mappings III models at SMC and LMC metallic- 
ity and assume that all the observed He n A4686 emission 
comes from shocks. We include emission from the shock's 
precursor, pre-shock gas that is photo ionized by UV an d 
X-ray emission produced in the shock ((Alien et alfe OOS) . 
We then derive the intrinsic photoionized line ratios us- 
ing 

A A -A s &X-fc (4) 



H Hq — H s 



X - 1 



where 



and 



X 



^ HeIIA4686 ^ 
^ HcIIA4686T 



(5) 



H 



HeIIA4686 = HcIIA4686 s 



(6) 

In the above equations, A is the emission line flux of in- 
terest, H is the H/3 line flux, the subscript o indicates 
observed values, and the subscript s indicates the mod- 
eled shock values. In addition to the two metallicities, we 
use the modeled shock values for magnetic field strengths 
of 0.5-10 fiG and velocities from 300-1000 kms" 1 . Fig- 
ures H2] and [TJ] compare the corrected photoionized line 
ratios for the two youngest Peas with the Starburst99 
and CLOUDY models described in f3~2~Tl 

The models show that when the additional flux con- 
tribution of shocks is taken into account, the optical 
depth implied by the nebular line ratios decreases. As 
the shock contribution to lower ionization lines, such as 
[On], increases, the intrinsic [O ill] /[Oil] ratio increases, 
and the Peas move into the density-bounded, optically 
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Fig. 12. — The panels show the stacked spectrum of the high [O m]/[0 II] ratio Peas for the spectral regions near shock lines. [Ne v] A3346, 
[Nev] A3426, [Fev] A4227, and [Mgl] A4571 are not detected. The positions of these wavelengths are indicated. [Ol] A6300 is shown in the 
bottom panel. 



thin regime. Furthermore, the presence of shocks can 
raise the Ha/H/ 3 ratio, leading to an o verestimate of the 
dust extinction (|Nakaiima et al.ll2012[) . If the dust abun- 
dance is overestimated, the Peas' [On]/H/3 ratios should 
be lower, and their [Oiii]/[Oii] ratios should be higher. 
This change would move the Peas farther into the opti- 
cally thin regime, in a direction opposite to the redden- 
ing arrow in Figures [6]|9] and IT3l[T4l The inferred optical 
depth of the Peas depends strongly on the inner nebu- 
lar radius and shock properties assumed, however. In 
particular, lower velocity shocks have a larger effect on 
the inferred optical depth. The shock emission also de- 
pends on the density, which may differ from the density 
of the photoionized gas. The SMC and LMC shock mod- 
els of I Allen ~a l. (2008) have a fixed density of 1 cm~ 3 . 
The Mappings III solar metallicity models indicate that 
higher density shocks have less effect on the Peas' ob- 
served line ratios. If the density of the shock-heated gas 
is greater than 1 cm -3 , the inferred optical depth would 



be greater than that implied by Figures Q2HH] An ad- 
ditional uncertainty comes from the modeled magnetic 
field strength, which is somewhat weaker than the Peas' 
magnetic fields. 

Although the presence of shocks shifts the Peas to- 
ward lower optical depths, the galaxies' line ratios still 
appear offset to higher [On]/H/3 ratios and/or lower 
[Oiii]/[Oii] ratios by factors of ^1.5-2 relative to the 
n < 1000 cm~ 3 model grids (see e.g., Figure IT51) . One 
possible explanation for this offset is that higher density 
models are more appropriate for the extreme Peas. If 
shocks are indeed present, they may contribute a sub- 
stantial fraction of the observed [S n] emission, which 
may lead to incorrect density estimates ($2j. The exis- 
tence of lower ionization nebulae within the Peas may 
also resolve the discrepancy. If low ionization nebu- 
lae contribute a fraction of the observed [O n] emission, 
they could account for the higher [On]/H/3 and lower 
[O in] /[On] ratios relative to the grids. Adopting lower 
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Fig. 13. — The grid is described in Figure [6] The figure shows the CLOUDY and Starburst99 grids for njj = 300 cm -3 and a filling 
factor of 0.1. The black cross indicates the observed ratios of Pea J101157. The la errors are indicated by the cross in the lower left of the 
plot. The squares and X's indicate the line ratios with the shock contribution subtracted. The squares correspond to shock models at LMC 
metallicity, and the X's indicate SMC metallicity. The symbol size indicates the magnetic field strength; large symbols correspond to 10 /iG 
models, intermediate-size symbols correspond to 3.23 fj,G, and small symbols correspond to 0.5 fiG. Symbol colors denote shock velocity: 
blue for 300 kms - 1 models, green for 500 kms" 1 models, and orange for 1000 kms -1 models. Accounting for a shock contribution to 
the Peas' emission decreases their inferred optical depth. 



reddening corrections or higher metallicities for the Peas 
could also lessen the offset; however, the uncertainties in 
these quantities are too low to completely explain the dis- 
agreement. Finally, as mentioned in £13.21 a population 
of hot, low- metallicity stars could provide an additional 
means of increasing the [O in] /[On] ratio. Given the 
current lack of evidence for such stars, however, we con- 
clude that the presence of multiple nebulae and higher 
densities in the star-forming regions are the more likely 
explanations. 

Despite these uncertainties, the CLOUDY models 
show that including shock emission lowers the inferred 
optical depths of the Peas, in some cases dramatically. If 
the shocks are due to SN activity, the Peas must be old 
enough for such events to take place. As a test, we con- 
sider Pea J 10 1157, the highest EW Pea and potentially 
the youngest. Our results for the other extreme Peas are 
comparable. From the observed Hell and the modeled 
Hell/Ha ratios, we estimate how many SNe are neces- 
sary to account for the He II emission. Assuming all the 
Hell is due to shocks, the modeled Hen/Ha ratios imply 
that shocks provide 5.6 x 10 40 -4.5 x 10 41 ergss -1 of Pea 
J101157's Ha luminosity. The Ha lumin osities of SNRs 
range from ~ 10 36 to ~ 10 38 ergss -1 (jde Griis et al.l 
2000). If we assume an intermediate SNR luminosity 
of 10 37 ergss -1 , then Pea J101157 would need 5600- 
45,000 SNe (—4-50% of its O stars) to reach the necessary 



Hell and Ha luminosities. While the former value could 
be reached in less than 4 Myr following an instantaneous 
burst, the latter value requires more than 6 Myr. How- 
ever, if the SNRs are extremely bright (10 38 ergss -1 ), as 
might be expected for a young burst, than the required 
number of SNe drops to only 560-4500, 0.4-5% of Pea 
J101157's O stars. The first SNe at -3.5-4 Myr could 
provide the necessary Ha and Hell luminosities. The 
ages for the continuous star formation case are similar; 
5000 SNe would go off within the first 5 Myr. 

One example of a Pea-like galaxy that may have 
shock- induced Hen emission is II Zw 40, a nearby 
BCD that shares many properties with the Green 
Peas. It has a metallicity of 0.2Z Q and an H/3 EW of 
2 72 A (|Guseva et al.l 120001) . suggesting a similar age to 
the Peas. In addition, II Zw 40 has detectable WR fea- 
tures, including the blue bump and Si in A4565, but the 
blue bump is dominated by the m uch stronger nebula r 
emission of Hen, [Fein], and [Anv] ()Guseva et al.l2000|) . 
The much weaker strength of the WR features compared 
to the nebular emission shows that these features could 
easily be below our detection limits. Based on the simi- 
larities betwee n He H-emitting gala xies with and without 
WR features, iGuseva et al.l (|2000l ) suggest that shocks 
generate some of the nebular He n emission, even in WR 
galaxies. The presence of weak [Fe v] A4227 emission and 
high- velocity emission-line wings in II Zw 40 supports the 
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Fig. 14. — The grid is described in Figure[6] and the symbols are described in Figure [T3l Both panels show the CLOUDY and Starburst99 
grids for n[j = 1000 cm" 3 . The model in the top panel uses a filling factor of 0.1, and the bottom panel uses a filling factor of 0.01. The 
black cross indicates the observed ratios of Pea J121903. 
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shock hypothesis (|Thuan fc Izotovl 120051 ). Although II 
Zw 40 is unlikely to be optically thin (e.g.. Uov &: Lesteil 
l!988tlBeck et al.ll2002l ). it shows that shock ionization of 
Hen is feasible at the Peas' metallicities and ages. Thus, 
even if WR stars are present in the Green Peas, shock 
heating may still be the primary origin of the high ion- 
ization nebular lines. 

The above analysis demonstrates that shock-heating 
from the earliest SNe could account for the observed 
Hen emission. In addition, shocks due to a galaxy 
merger or gas infall could occur even earlier, and younger 
ages for the Peas could be possible. If shocks are present, 
the Peas' high [Oiii]/[Oii] ratios may indicate a low op- 
tical depth. 

4. CONSEQUENCES FOR THE OPTICAL DEPTH 

The source of the high ionization emission in the ex- 
treme Green Peas could be WN9 stars, early-type WR 
stars, chemically homogeneous O stars, shocks, or some 
combination of these sources. Each of these sources has 
different implications for the escape of ionizing radiation, 
which we now discuss. 

If the He n emission is a stellar feature or due to stel- 
lar photoionization, the Peas' optical depth is ambigu- 
ous. Low density models, with n < 100 cm~ 3 indicate a 
high optical depth. At higher densities, consistent with 
the Peas' [Sn] ratios, whether or not the Peas are opti- 
cally thin depends on the inner nebular radius and the 
filling factor. For instance, decreasing the filling factor 
of the dense gas decreases the inferred optical depth (see 
Figures [6] and El). 

Increasing the inner radius lowers the ionization pa- 
rameter; the high [O in] /[On] ratios we observe would 
then have to be due to a low optical depth. Unlike 
the CLOUDY models, the emission from the Peas is not 
likely to originate from a single H II region. However, 
if the emission is dominated by one cluster or several 
similar clusters, we can estimate which CLOUDY model 
geometries may be most reasonable. 

The stellar winds of WR st ars may bl ow cavities of 50 
pc or more m the ISM fe.g.. lQg^ll996h . and SN explo- 
sions will expand the superbubble cavity further. The 
broad wings in the Peas' Ha emission may already indi- 
cate that stellar feedback is evacuating g as around the 
star cl us ters and driving gala c tic w inds (jAmorm et al.l 
I2012H ). IStrickland fc Stevens! (|1999| ) simulate a super- 
bubble around a 10 5 M Q cluster at 0.25 Z , expand- 
ing into an ambient medium with density n = 100 
cm~ 3 . The superbubble cavity reaches a radius of 50 
pc within 3.5 Myr and expands to 70 pc within 4.5 Myr. 
The radius scales with cluster mass as R oc (M /n) 1 / 5 
(jCastor et alj 119751: IStrickland fc Stevens! H999I ). so for 
our 10 6 M Q Starburst99 simulation, the cavity radii 
should be 79 and 111 pc, respectively. If we increase 
the density by a factor of 10, as appropriate for Pea 
J121903, the radii remain 50-70 pc. Therefore, for our 
10 6 M Q starburst CLOUDY models, the models with in- 
ner radii greater than 50 pc may have the most realistic 
geometries. Models with these inner radii and n > 100 
cm~ 3 generally indicate a low optical depth (see Fig- 
ure©. 

The exact optical depth of the Peas also depends on 
their age and the resulting ionizing spectrum. For in- 
stance, Figure [6] shows that for a given inner radius, a 5 



Myr-old model gives a slightly lower optical depth than 
a 4 Myr-old model. The ionizing spectrum of chemically 
homogeneous O stars is highly uncertain, and the exis- 
tence of 60,000 K homogeneous stars has not been obser- 
vationally confirmed. We cannot determine the optical 
depth of the Peas if their ionization is due to a chem- 
ically homogeneous population. Thus, while we cannot 
assess the precise optical depth of the extreme Peas, a 
low optical depth remains plausible, particularly for the 
Peas with the highest densities. 

As discussed in M3.51 including a shock contribution to 
the Peas' line fluxes decreases estimates of their optical 
depth. This same situation is appare nt in the s pectra 
of DEM L301, an LMC superbubble. IQev et al.l pOOOl ) 
show that a combination of shock-heating and density- 
bounding explain DEM L301's emission, although at 
first glance , the line rat i os app ear to indicate radiation- 
bounding. iVoges et all (|2008[) confirm that DEM L301 
is indeed optically thin. 

Not only do SNe affect the nebular emission, but they 
may actually promote a low optical depth. As SN-driven 
bubbles expand in the ISM and eventually blow out of 
the galaxy, th ey create low-densit y passages that aid 
photon escape ([Clarke fc Oe v 2002). Like Lyman-break 
galaxies, the Peas have large SFRs and spatially dis- 
trib uted star formation, which enhance the ISM poros- 
ity (|Clarke fe; Oevll2~002f h In addition, with their young 
ages, the Peas may be at the optimal time for the escape 
of Lyman-continuum radiation; they are old enough for 
SNe and stellar winds to begin to reshape the ISM, but 
young enough to possess large numbers of UV-luminous 
massive O or WR stars. 

The CLOUDY models for the WR and shock-heating 
scenarios at n — 100 cm~ 3 are generally optically thick 
or borderline optically thin. Models at higher densities 
favor lower optical depths but could still be consistent 
with optically thin nebulae if the inner radii are suffi- 
ciently small and shocks are negligible. Nevertheless, at 
all densities, the Peas may be more optically thin than 
the models suggest. For instance, optically thick clumps 
tend to dilute the signature of density -bounding (e.g., 
IQev et al.l[2000l : lGiammanco et al.ll2005l ). The CLOUDY 
models further assume that the emission originates from 
a single nebular region or identical nebulae. In reality, 
the emission is the sum of the emission from multiple 
nebulae, which may have different levels of ionization. 
Low-ionization emission from other nebulae or a com- 
panion galaxy could lower the observed [O in] /[On] ra- 
tios and cause the gas to appear more optically thick 
than it actually is. 

The extreme Peas may therefore have high escape frac- 
tions for ionizing radiation. However, further observa- 
tions are needed to diagnose their optical depths and 
distinguish among the potential ionizing sources. We 
are obtaining deeper spectra of the extreme Peas to con- 
firm or refute the existence of shocks and to constrain 
the WR population. Observations of massive stars in 
low-metallicity galaxies are also critical to understand 
homogeneous stellar evolution and the stellar spectra of 
low-metallicity O and WR stars. 

5. SUMMARY 

The Green Peas are a rare class of compact, emission- 
line galaxies at z = 0.1 — 0.3 with unusually powerful 
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[O in] emission and the potential for a high escape frac- 
tion. Through their high ionization emission, extreme 
star-formation intensities, and clumpy morphologies, the 
Green Peas provide a glimpse of massive star feedback 
and ionizing photon propagation in high rcdshift galax- 
ies. While the large [Oiii]/[Oii] ratios of the Peas may 
indicate a low optical depth, high ionization lines such 
as He II and [Ar iv] in the SDSS spectra show evidence 
for strong, high-energy emission. Using stellar popula- 
tion and nebular photoionization models, we have inves- 
tigated the ionizing sources and optical depths for six of 
the highest ionization Peas. 

The six extreme Peas are young and powerful star- 
bursts, with high electron temperatures of ~ 15, 000 K 
and densities of 100-1000 cm~ 3 . Their high Balmer- 
line EWs indicate young ages of 3.7-5.1 Myr, according 
to Starburst99 instantaneous burst models. We detect 
the He I A3819 line in emission in several objects and in 
stacked spectra, which corroborates the young ages. The 
A3819 line sets an upper limit of 3 Myr on the ages for 
an instant aneous burst and 10 Myr for continuous star 
formation (jGonzalez Delgado et al.|[T999[ ). 

The Peas' young ages and the presence of He n A4686 
emission set strong constraints on the possible ioniz- 
ing sources in the Peas. No WR features are visible in 
stacked spectra of the Peas. However, we tentatively de- 
tect WR features in Peas J101157, J145735, and J030321. 
The narrow FWHMs (3-5 A) of the He II emission suggest 
a nebular origin or a narrow WNL st ellar feature. 

Usin g th e WNL line luminos ities of iBrinchmann et al.l 
(2008J) and lGuseva etal] (pOOOT) , we calculate the number 
of WNL stars necessary to produce the observed A4686 
emission as a stellar feature. The high resulting WR/O 
ratios and the 3er Si in A4565 detection limits rule out the 
required WNL population for the Z < 0.2 Z Q case, but 
a WNL population with Z > 0.2 Z Q is possible. How- 
ever, a higher-metallicity WNL population may not be 
consistent with the Peas' low measured metallicities. 

Alternatively, the observed He II may result from pho- 
toioniza tion by hot stars. CLOUDY photoionization 
models ([Ferland et al.l 119981 ) of single-star H n regions 
demonstrate that ordinary O stars are not hot enough 
to provide the necessary He + -ionizing flux. CLOUDY 
models using early-type WR stars, with core tempera- 
tures above 70,000 K, fit most of the observed l ine ra tios. 
From the ionizing photon rates of lSmith et al.l ()2002D . we 
estimate the number of WNE and WCE stars needed to 
produce the observed Hen. The resulting WR/O ra- 
tios are consistent with Starburst99 instantaneous burst 
models, but the lack of broad He n emission in most ob- 
jects eliminates the WCE stars as the dominant cause 
of the emission. Starburst99 and CLOUDY models like- 
wise show that only 4-5 Myr starbursts, which contain 
a substantial WR population, can simultaneously match 
the Peas' [Om]/[Oii] and Hen emission. 

Putative chemically homogeneous O stars could ex- 
plain Hen emission in galax ies without WR features 



plain rle II emission m galax ies without WK features 
(jShirazi fc Brinchmannl 20121) . These stars may form 
at low metallicity (e.g., (Maederl Il987l: lYoon et all 120061) 
and can theoretically reach effectiv e temperatu r es muc h 
higher than ordinary O stars ()Brott et al.l l2011al ). 
CLOUDY single-st ar H II region mod els with the 60,000 
K stellar model of IKudritzkl (|2002D demonstrate that 



such hot O stars could account for the observed Hen. 
In addition, the presence of extremely hot O stars 
may be necessary to match the observed [On]/H/3 and 
[O in] /[On] ratios of the Peas. More research is needed 
to determine whether homogeneous O stars exist, how- 
ever, and how many ionizing photons they produce. 

We rule out an AGN contribution to the Hen emis- 
sion at the 10% level, using the diagnostic diagram of 
Sh irazi &: Brinchmannl (j2012). Likewise, we find that 
HMXBs cannot be t he source of the H en ionization. 
Using the models of ILinden et al.l (|2010), we calculate 
the expected number of luminous HMXBs in the Peas. 
Even assuming that each HMXB is a ULX, with the 
high nebular Hell A 4686 luminosity of Holmberg II X-l 
(|Kaaret et al.|[2004l ) , the resulting emission is an order of 
magnitude too low. 

Fast, radiative shocks are an alternative mecha- 
nism for the product ion of He II (|Guseva et al.l 120001 : 
iThuan fc Izotovir2005D . and the existence of fast shocks 
is supported by the high velocity Balmer-line wings 
(jAmorin et al.ll2012b[) and [O i] and [N i] emission in the 
Peas. To account for the contribution of shocks to the 
observed emission, we use the line emission predict ions 
of the Mappings III shock models ()Allen et al.ll2008l ) for 
different magnetic field strengths and velocities. Assum- 
ing shocks generate all of the Hen emission, we sub- 
tract the shock contribution from the observed line ra- 
tios. After correcting for shock emission, the nebular 
line ratios indicate a lower optical depth, particularly 
for slower shocks. Based on the Ha luminosities from 
the shock models, the first SNe, exploding 3.5-4 Myr af- 
ter an instantaneous burst, could produce the observed 
He II emission. Younger ages could be possible if infalling 
gas or mergers are the cause of the shocks. 

Although the precise origin of the He n emission is still 
uncertain, the CLOUDY models show that the Peas may 
be optically thin. Models with higher densities and larger 
inner radii imply lower optical depths. If the Hen is due 
to WR stars, models with inner radii > 50 pc and den- 
sities ~ 100 cm -3 support both optically thick and op- 
tically thin scenarios, while models with the same inner 
radii and n = 1000 cm -3 support a low optical depth. 
Including the effects of shocks lowers the estimated op- 
tical depth further, and the presence of density inhomo- 
geneities or emission from multiple nebulae may also lead 
us to overestimate the optical depth. 

The low levels of extinction, intense star formation, 
and young ages of the Green Peas make them ideal candi- 
dates for escaping ionizing radiation. The nebular emis- 
sion of the Peas with the highest [O m]/[0 n] ratios, their 
high densities, and the likely presence of shocks in these 
objects suggest that they may be leaking ionizing pho- 
tons into the IGM, and their ages may optimize photon 
escape. In these 3-4 Myr starbursts, stellar feedback is 
likely beginning to carve out passages in the ISM for the 
escape of ionizing photons, while large numbers of mas- 
sive stars are still available to produce Lyman-continuum 
photons. However, the data are still inconclusive, and 
future observations are necessary to establish the escape 
fraction of the Green Peas. These galaxies offer a rare 
opportunity to study high-redshift star formation condi- 
tions, the effects of stellar feedback on the ISM, and the 
propagation of ionizing photons from star clusters to the 
IGM. 
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